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We examine the conditions of validity for the Lifshitz-Matsubara sum formula for the Casimir 
pressure between magnetic metallic plane mirrors. As in the previously studied case of non-magnetic 
materials (Guerout et al, Phys. Rev. E 90 042125), we recover the usual expression for the lossy 
model of optical response, but not for the lossless plasma model. We also show that the modes 
associated with the Foucault currents play a crucial role in the limit of vanishing losses, in contrast 
to expectations. 
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I. INTRODUCTION 


The comparison of experimental measurements of 
Casimir force with theoretical predictions remains a mat¬ 
ter of debate HH3|. For experiments performed with mir¬ 
rors covered by thick layers of gold, optical properties of 
the metallic mirrors should be deduced from tabulated 
optical data II 0 and extrapolated to low frequencies 
by using the Drude model, with a finite ohmic dissipa¬ 
tion rate 7. However experimental data EUZ] appear to 
be in better agreement with the so-called plasma model 
which assumes 7 to vanish, in clear contradiction with 
the well established fact of a finite static conductivity of 
gold. This “Casimir puzzle” remains to be solved. 

In a previous paper [ 5 ] , we re-examined the conditions 
of validity of the Lifshitz formulas used to calculate the 
Casimir pressure. We studied in a very cautious way 
the positions of the modes of the system in the complex 
frequency plane, and identified a previously unsuspected 
problem in the use of the Lifshitz-Matsubara sum formula 
for the lossless plasma model. As it is usually applied, 
this formula neglects the contribution of the modes asso¬ 
ciated with Foucault currents, although the latter tends 
to a non vanishing limiting value when 7—^0. While 
this work did not solve the aforementioned Casimir puz¬ 
zle, it cured the discontinuity of the calculated Casimir 
pressure at 7 = 0 . It was also interesting from a pure 
theoretical point of view, as it shed new light on the sub¬ 
tleties of the application of Cauchy’s residue theorem in 
the context of the calculation of the Casimir pressure. 

The Casimir puzzle is not only apparent in experiments 
performed with gold. It has also been confirmed in recent 
experiments involving magnetic materials [TO) 111 ] where, 
again, experimental measurements appear to agree with 
predictions based on the non physical lossless plasma 
model, rather than the much better motivated lossy 
Drude model. This situation has pushed us to extend 
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our study to include the case of magnetic metallic mir¬ 
rors. This requires that the mirror’s optical properties 
be also described by a frequency-dependent permeabil¬ 
ity r(oj), besides the more commonly studied frequency- 
dependent permittivity e(w). This situation leads to a 
whole new structure of low frequency modes with a richer 
phenomenology than in the non-magnetic case. 


II. OUTLINE 

The Casimir pressure between plane mirrors can be 
calculated equivalently as a Lifshitz integral over real 
frequencies oj or as a Lifshitz-Matsubara sum over imag¬ 
inary frequencies *£ n . The two formulations are mathe¬ 
matically related via the application of Cauchy’s residue 
theorem, and the precise conditions of validity of this 
equivalence are discussed in [8]. The optical properties 
of the mirrors are described by reflection amplitudes for 
the two polarizations TM and TE, written in terms of 
permittivity and permeability functions by using Fresnel 
laws. 

We begin with non-magnetic materials, for which the 
reduced permittivity function is written as 

e(w) = 1 + x(w) , X'M = x'M + ix" M , (1) 

where x'( w ) an d x!'( u ) are real f° r a ^ real frequencies 
oj. As the dielectric response is causal, the function x( w ) 
has an analytic continuation x( z ) i n the upper half-plane 
which decays fast enough, at least in l/|z|, so that inte¬ 
gration contours can be closed at infinity. It follows, for 
smooth enough functions, that they obey the Kramers- 
Kronig relations 


X{z) 



x\{x) 


d;r. 


( 2 ) 


for z a complex number in the upper half-plane. More 
generally, causal response functions obey dispersion re¬ 
lations which modify these Kramers-Kronig relations 
(examples discussed below). 




2 


As the response function is real in the time domain, 
x(z) obeys the Schwartz reflection principle along the 
imaginary axis x*(~) = x{~ z *)- For z = lying on the 
imaginary axis, one gets the familiar relation 
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The knowledge of the dissipative part xM w ) = e"M 
of the permittivity is sufficient for constructing x(*£) = 
e(*£) — 1 . The latter is deduced from the optical data 
tabulated over some interval of frequencies from far in¬ 
frared to ultraviolet in the best scenario, and extrapo¬ 
lated to low and high frequencies HUS]- In the following, 
we focus the discussion on the contribution of conduc¬ 
tion electrons, dominant at low frequencies. Of course, 
the contribution of bound electrons has to be included as 
well in the full analysis. 

For metals, the experimental data for xMM have to 
be extrapolated at low frequencies by using for example 
the Drude model 
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<t 7 M is the conductivity, w p the plasma frequency and 7 
the dissipation rate which determines ohmic losses [12] . 
The Drude model is the simplest one to match the im¬ 
portant fact that metals have a finite static conductivity 

w 2 

<t 7 (u; —> 0 ) = — ( 5 ) 

7 

and thus to be compatible with Ohm’s law. We note 
that the function x 7 obeys the relation ([3]) used in many 
papers as a starting point of the evaluation of Casimir 
forces. The function xx"{x) appearing in the numerator 
in the integral © is the real part cr'(x) of the conduc¬ 
tivity and it is regular everywhere. 

The lossless plasma model is obtained by setting 7 = 0 
in the Drude model 

XoM = —— , <Jq (w) = (6) 

— UxJ —VjJ 


This simplified model can be a fair approximation of the 
contribution of conduction electrons at large frequencies 
u) 7. However it contains no dissipative part, thus 
failing to reproduce optical data at low frequencies, and 
it also misses the fact that metals such as gold have a 
finite static conductivity ([ 5 ]). The plasma model ^ does 
not obey the relation ®, as is obvious from the mere 
fact that xXq{x) = a q(x) vanishes, so that ^ implies 
Xo(*0 = 0, in contradiction with ([6]). 

In more mathematical terms, this defect can be at¬ 
tributed to the fact that xo has a double pole at the 
origin or equivalently, that 00 has a simple pole at the 
origin. A proper application of Cauchy’s theorem then 
leads to the following expression 


Xo(z) 
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where the integral contribution is in fact zero as Xo = 0. 
In other words, Kramers-Kronig relations do not have 
any physical content for models which do not include 
dissipation. For the lossless plasma model, dispersion 
relations are written with subtractions [ 5 ] , which are just 
modifications of Kramers-Kronig relations. 

In our minds, the agreement of measurements with pre¬ 
dictions of the plasma model constitutes a problem to be 
solved, the Casimir puzzle mentioned above. The lossless 
plasma model, with susceptibility xoMi does not match 
the optical and electrical properties of gold, and it can¬ 
not describe correctly the Casimir pressure between two 
metallic plates. As already stated, the plasma model can 
only be considered as an effective model at high frequen¬ 
cies u 7. Considered in this manner, it has to be 
defined as the limit of the Drude model when 7 —>■ 0 

X (w) = a ( w ) = Up . , 77 —S> 0+ . (8) 

— VjJ T] — VjJ 

This definition takes care of the problem shown by Xo( ta ) 
at ui —> 0 by a proper specification of the function there, 
in the spirit of causality arguments. It gives a suscepti¬ 
bility Xr;( w ) which differs from x 0 ( c * ; )i with the difference 
identified as what is left of dissipation in the vicinity of 
u> = 0 when 77 —► 0 + . A careful use of the theory of 
distributions [ 13 ] leads to an expression of this difference 

w 2 

X77M = XoM - otcu^'M > XoM = -J 2 • (9) 

The discussions in [8] prove that predictions based on Xv 
are then obtained by continuity of those of the Drude 
model x 7 - In contrast, the predictions based on the loss¬ 
less plasma model xo, as advocated for example in fl 4 l . 
show a discontinuity with those of the Drude model that 
we will interpret as reflecting the difference ©>• 

Let us now focus on the extension of the results in [ 5 ] 
to situations involving magnetic mirrors. We calculate 
the pressure between two metallic parallel plane mirrors. 
The mirrors are finite-thickness slabs but we consider the 
thickness to be large enough so that the mirror’s reflec¬ 
tion coefficients are indistinguishable from Fresnel ampli¬ 
tudes. The dielectric properties of mirrors are described 
by a Drude model at low frequencies. We model the 
magnetic response by a frequency-dependent permeabil¬ 
ity as in [IT] 


/r(w) — 1 + 


MO) -1 

1 - iuj/LO m 


( 10 ) 


This relaxation form is appropriate for describing the 
spin rotational component of the magnetization. /r(w) 
monotonically decreases from /i(0) to 1 over a character¬ 
istic frequency w m and it has a simple pole at z = —ioJ m 
in the complex plane. 

The characteristic frequency is by far the lowest 
frequency in the problem. In El it is argued that ui m lies 
in a frequency range of the order of 10 5 - 10 9 Hz. We may 
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emphasize at this point that this is a further complication 
for the plasma model scenario. For non magnetic mirrors, 
the parameter 7 is indeed the lowest frequency in the 
problem. In contrast, for magnetic materials its physical 
value is much larger than that of u> m and it makes even 
less sense to consider 7 —> 0. In the following, we will 
consider different values for 7 , but always stress that its 
physical value obeys 7 w m . 

In order to specify the other parameters, we will use 
particular values in the models of permittivity and per¬ 
meability which match the non-magnetic and magnetic 
metals used in the experiments, that is to say gold (Au) 
and nickel (Ni). For gold, we will use Huj p = 9 eV, 
hjAu = 35 meV and /r(0) = 1; For nickel, tuo p = 4.89 
eV, = 43.6 meV, /x(0) = 110 and hoj m ~ 0.1 neV. 

The magnetic permeability enters the definition of the 
Fresnel amplitudes as 


„TE 
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pk z - K z 
pk z + K z 


TM 
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ek z - K z 
ek z + K z 


K z = 



( 11 ) 


In the following we focus on the TE polarization. The 
Fresnel amplitude r™ when extended to complex fre¬ 
quencies z = ui + i£ vanishes as z at the origin for non¬ 
magnetic materials. For magnetic materials, it behaves 
as z° at the origin instead. We use notations inspired 
from our previous work [5] and write the Casimir pres¬ 
sure P c calculated using the Lifshitz formula as 


p c 

/£ M 



d uj 
27r 


Re 14]. 4 


hk z f^)c M, (12) 


(rl(uj)) 2 e 2lk * L 
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C(uj) = coth 
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where L is the distance between the mirrors, T the tem¬ 
perature and c the polarization. We also recall the 
Casimir pressure P CM calculated using the Lifshitz- 
Matsubara sum formula as 


P CM = —2 k B T (is) 

k,s n 



where the *£„ are the Matsubara frequencies. In this 
context of the the different models of permittivity used 
to describe metallic mirrors, we have shown in [ 8 j that 
the application of Cauchy’s residue theorem leads to 

P~-Ef (14) 

= lim / 

£_>0+ Jr\[—e,e] 



so that P CM coincide with P c except if the spectral den¬ 
sity Re[p^(w)] contains a Dirac S(u>) at the origin. Inde¬ 
pendently of which Lifshitz formula is used, one can use 


the different susceptibilities introduced before. The sus¬ 
ceptibility Xv( u ) bring a distribution leads to an unde¬ 
fined quantity. On the other hand, we define the Casimir 
pressures Pq and P 7 as the results of the formulas ( 12 ) 
and (13) when using the susceptibility Xo(w) and x 7 (uJ). 
Those are well defined quantities and one can study the 
limit 


P v = lim P-y. (15) 

7—>-0 

In [5], we have shown that this limit exists and that 

Pf=P 0 C - P£ M ^P£ M . (16) 

Unless otherwise stated, we will work in the following 
with the quantity P 7 for which the two Lifshitz formulas 
coincide. 

We now recall in more detail the main results of our 
previous work (5j which uses an elegant extension of the 
argument principle. In a system of two non-magnetic 
lossy mirrors the function p™(z) behaves as z 2 z~ 1 = z 
at the origin. Therefore, it makes no contribution to the 
Lifshitz-Matsubara sum formula for the Matsubara fre¬ 
quency £ = 0 (the Lifshitz-Matsubara sum formula picks 
up a non-vanishing contribution when 4 (~) behaves as 
(z — £„) _1 at a Matsubara frequency £„). At the same 
time, the modes associated with the Foucault currents 
lie in an interval £ £ [— 7 , — 7 ], 7 > 7 > 0, on the neg¬ 
ative imaginary axis of the complex plane. The point 
z = — *7 is an accumulation point. Thus, there is an in¬ 
finite number of poles (modes) and zeros in the interval 
£ £ [- 7 ,- 7 ]. Nevertheless, one can evaluate the quan¬ 
tity 

N = £d z lo gP r(z)^ (17) 

on a positively-oriented closed contour C surrounding the 
interval £ £ [— 7 , — 7 ] to find N = —2. This means that 
Pk E (*) possesses two more poles than there are zeros in 
this interval. As 7 —> 0, the interval [— 7 , — 7 ] collapses 
on the simple zero at the origin leading finally to the 
function p r £ E (z) behaving as z~ x there. The Lifshitz- 
Matsubara sum formula now picks up a non-zero contri¬ 
bution from the TE polarization and the Casimir pressure 
is therefore discontinuous. But as 7 —> 0 the collapse of 
the Foucault modes in the interval [— 7 , — 7 ] gives rise to 
a Dirac S(co) contribution in the function Re[p™( w )]- 

As mentioned before, this contribution is omitted in 
the Lifshitz-Matsubara sum formula and must be taken 
into account explicitly. When this is done it is seen that 
the contribution from the S(co) at the origin exactly can¬ 
cels out the Matsubara pole contribution there and the 
Casimir pressure is therefore continuous. We then con¬ 
clude that the plasma prescription applied to the Lifshitz- 
Matsubara sum formula fails in taking into account the 
contribution from the Foucault modes which persists 
even in the limit of lossless metals. For non-magnetic mir¬ 
rors, the transverse wavevector k is a spectator through¬ 
out this process only influencing the value of 7 on the 
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negative imaginary axis. In particular, we have shown 
that the contribution from the Foucault modes to the 
Casimir pressure is always repulsive for all k for non¬ 
magnetic mirrors. 

In the following, we study the position of the low fre¬ 
quency modes of the system containing magnetic mirrors 
as 7 —> 0 . For mixed magnetic-non-magnetic situations 
both relaxation rates, e.g. 7a u and 7Ni, go to zero at the 
same time. An asymptotic regime is achieved when 7 is 
smaller than the smallest frequency in the problem. In 
particular, we consider in the following 7 < w m which is 
not in the physical region. 


III. THE AU-NI SYSTEM 


We have mentioned that for the Drude permittivity, 
the point z = —*7 is a pole and an accumulation point. 
Similarly, for the permeability given by eq. ( 10 ) the point 
z = —iLOm is also a pole and an accumulation point. The 
permeability n(z), extended to the complex plane, gives 
rise to magnetic modes which lie on the negative imag¬ 
inary axes in an interval £ € [—w m ,— w m ], uj m > ui m . 
Both points z = —ioj m and z = —i'y are zeros of K z . In 
the limit 7 /uj m < 1 we have 


^ k 2 c 2 + Wp/i(0) 

k 2 c2+LJ 2 ’ 

k 2 c 2 

^ ^ k 2 c 2 + Wp/Lt(0) 


( 18 ) 

( 19 ) 


The Fresnel TE amplitude given by eq. © now pos¬ 
sesses two zeros on the imaginary axes, which we de¬ 
note by {£o~i£(^}' The position of those zeros depend 
on all the parameters of our system. Notably, they de¬ 
pend on k = |k|. Nevertheless, they satisfy some general 
properties: one of these zeros always lies on the positive 
imaginary axis > 0, the other zero always lies on the 
negative imaginary axis and we have —< —7^] 
Finally, the function Pk E (~) behaves as zz°z ~ 1 = z° at 
the origin so that the TE polarization at £ = 0 does not 
contribute to the Casimir pressure. We show in fig. [l] 
the schematic representation of the low frequency inodes 
of the Au-Ni system, as poles and zeros of the function 
p™(z), in the vicinity of the origin of the complex fre¬ 
quency plane. The modes associated with the Foucault 
currents are, strictly speaking, both the Foucault modes 
from the gold and nickel mirrors. As such, the value 7 in 
the figure is really max( 7 A u , 7 Ni) and 7 = min( 7 A u , 7 Ni)- 
This whole combined structure has N = —2 as mentioned 
before. A calculation of N on a contour which encloses 
the interval £ £ [—Cj m , — uj m ] gives IV = —1 meaning that 


1 r 'k = 1 and r™(—* 7 ) = —1. Since E R for 

£ E [— uim, — 7 ] this interval must contain at least one zero. 



FIG. 1. Schematic representation of the poles and zeros of the 
function p™ in the vicinity of the origin for the Au-Ni system 
(mixed magnetic-non-magnetic mirrors). Zeros are indicated 
by crosses. Dashed intervals represent an infinite number or 
poles and zeros. Here we are considering 7 < u> m in order to 
study the limit 7 —> 0 . 


there is one more pole than there are zeros in this inter¬ 
val. 

From our previous work (Sj we know that the interval 
£ £ [—7, —7] collapses to the origin as 7 —>■ 0 . When 
dealing with magnetic materials, it turns out that one of 
the two zeros {£o";£o~} a ls° collapses to the origin. This 
leads at the limit 7 = 0 to a function p^ E (z) behaving as 
z" 1 at the origin, as it should. Which of the two zeros 
U 0 -,£ 0 + } collapses to the origin depends on the wavevec- 
tor k: it is easy to show that there is a remarkable value 


k 0 c 


yV(0Vp 

V^ 2 (0) - 1 


( 20 ) 


which corresponds to the condition r™( 7 = 0, w ^0) = 
0 . Depending on the value of k with respect to fco, one 
of the two zeros {£o"j£o^ } collapses to the origin as 

^7^o~l-fc2/fc2 (21) 

while the other zero tends to a finite value as 7 —>■ 0 . Let 
be the interval £ £ [—7, —7] supplemented by either 
^ or £+: 

: -7<e<-7U( = e 0 ± . (22) 

The collapse of either set S~~ or <S + to the origin is accom¬ 
panied in the function Re[pk E ( w )] by the appearance of 
a Dirac S(uj) distribution representing the non-vanishing 
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contribution of the Foucault currents at 7 = 0 for this 
particular wavevector k. Interestingly, when S~ or S + re¬ 
spectively, collapses at the origin the function Re[p™(w)] 
tends towards a positive or negative, respectively, Dirac 
6 (ui) at the origin. Therefore, when dealing with mag¬ 
netic materials the total contribution from the Foucault 
currents at 7 = 0, integrated over all fc, can in princi¬ 
ple be either repulsive or attractive. This is in contrast 
to what occurs with non-magnetic materials where this 
contribution was always repulsive (8]. 


Fresnel amplitude £(j“ and £(/ are now double zeros. The 
set of purely magnetic modes which lie in the interval 
£ £ [—w m , — w m ] now correspond to the modes of the two 
mirrors. As such, a calculation of TV on a contour enclos¬ 
ing this interval now leads to N = — 2 . The combined 
sets S~ and S + have N = 0 now. When they collapse at 
the origin, the function p™(z) still behaves as 2 -1 . 


V. TOTAL FOUCAULT MODES 
CONTRIBUTION 


IV. THE NI-NI SYSTEM 

The system of two identical magnetic materials inherits 
many of its analytic properties from the previous Au-Ni 
system with some slight differences however. We show in 


We show in fig. [ 3 ] the contribution of the Foucault 
modes as a function of the wavevector k for the three 
systems Au-Au, Au-Ni and Ni-Ni. Positive or negative 
parts, respectively, correspond to repulsive or attractive 
contributions. The total contribution to the Casimir 



FIG. 2. Schematic representation of the poles and zeros of the 
function p™ in the vicinity of the origin for the Ni-Ni system 
(identical magnetic mirrors). Zeros are indicated by crosses 
together with their multiplicities, poles by a dot. Dashed 
intervals represent an infinite number or poles and zeros. 

the schematic representation of the low frequency 
modes of the Ni-Ni system, as poles and zeros of the func¬ 
tion pJ E (x), in the vicinity of the origin of the complex 
frequency plane. First of all, the function p™(~) now 
behaves as z°z°z ~ 1 = z~ x at the origin so that the TE 
polarization now contributes to the Lifshitz-Matsubara 
sum formula at £ = 0 . This is represented by a simple 
pole which sits at the origin in fig. [ 2 ] This simple pole 
is handled in the usual way by choosing a contour of in¬ 
tegration which avoids it (see for instance fig. 1 for the 
TM polarization of our ref. 0). This leads to the usual 
contribution to the Casimir pressure at zero frequency 
which is weighted by a factor 1 / 2 . The zeros of the TE 
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FIG. 3. (Color online) Contribution to the Casimir pressure 
from the Foucault modes as a function of the wavevector k. 
From top to bottom, for the Au-Au, Au-Ni and Ni-Ni systems. 
Positive (red) parts correspond to repulsive contribution and 
negative (green) parts to attractive ones. 


pressure is the integral over k. This contribution is easily 
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calculated by setting the parameter 7 much smaller than 
all other frequencies in our problem so that the modes 
associated with Foucault currents, sitting at frequencies 
ui ~ 7, are well separated from all other modes. Then 
the quantity 

Jp g2Re[pT] (23) 

gives the contribution from those modes for a particular 
wavevector k. 

As mentioned before, in the case of purely non¬ 
magnetic materials (as exemplified by the Au-Au system) 
the contribution from the Foucault modes are always re¬ 
pulsive for all k. On the contrary, in the case of purely 
magnetic materials the contribution from the Foucault 
modes is almost always attractive (except for a negligi¬ 
bly small repulsive contribution at low k ). In the mixed 
Au-Ni system, the overall contribution from the Foucault 
modes is smaller in magnitude. In addition to that, there 
is now clearly a competition between repulsive and at¬ 
tractive contributions. In the example shown in fig. [3] 
the total contribution is slightly attractive. 

We present in fig. [4] a calculation of the Casimir pres¬ 
sure for Au-Au, Au-Ni and Ni-Ni plane mirrors separated 
by a distance L = 300 nm as a function of the relaxation 
parameter 7 used in the low frequency part of the Drude 
permittivity e 7 (i£). At 7 = 0, the triangles represent 


Pressure (mPa) 



FIG. 4. (Color online) Calculated Casimir pressures P 7 at 
T = 300 K for the Au-Au, Au-Ni and Ni-Ni systems of plane 
mirrors separated by a distance L = 300 nm (from bottom to 
top) as a function of the relaxation rate 7 . The squares are 
the calculations using relaxation rates given in the text. The 
triangles are Pq M using the Lifshitz-Matsubara sum formula. 
The circles are Pq using the Lifshitz formula which naturally 
take into account the Foucault currents. 

Pq M obtained using the Lifshitz-Matsubara sum formula 
whereas the circles represent Pq obtained using the Lif¬ 
shitz formula. The Foucault modes contribution shown 
in fig. 0 are directly apparent as the difference between 
the triangles and the circles. 

In the mixed case Au-Ni, it was noted in fTO] that 
no notable differences was seen in the experimental data 
and in the calculation between the Drude and the plasma 


prescription. This fact originates in a coincidental near¬ 
degeneracy between the two prescriptions in this case. 
The puzzle mentioned from the introduction is the fact 
that experimental data are in very good agreement with 
the triangles in fig. [4] even though those correspond to 
inconsistent calculations as we have shown. 

In addition to the direct comparison with experiments 
which measure the Casimir pressure, a recent work j!5] 
has shown that whether the calculation is performed us¬ 
ing the Drude or plasma models has an impact on the 
stability of actuation dynamics of microelectromechani¬ 
cal systems. For this reason, we show in fig. [5] the con- 

Foucault 
contribution (Pa) 



FIG. 5. (Color online) Contribution from the Foucault modes 
to the total Casimir pressure as a function of the distance L 
between the plates. From top to bottom for the Au-Au, Au- 
Ni and Ni-Ni systems. The inset shows the same data, in 
absolute value, on a log-log scale. 

tribution from the Foucault modes as a function of the 
distance L between the plane mirrors. For the Au-Au 
system, this contribution stays repulsive at all distances. 
However, for the Au-Ni and Ni-Ni systems, the contri¬ 
bution from Foucault modes stays attractive at all dis¬ 
tances. We note that for the mixed Au-Ni system, the 
Foucault mode contribution which was relatively small 
at L = 300 nm (see fig. [4]) becomes larger in magnitude 
than the Au-Au system at shorter distances (see the inset 
of fig. [5] which shows the contributions, in magnitude, on 
a log-log scale for the three systems). 

VI. CONCLUSION 

In this paper, we have confirmed that for both mag¬ 
netic and non-magnetic metals the contribution from the 
Foucault modes to the Casimir pressure between plane 
mirrors reduces to a non-vanishing Dirac 5(to) at vanish¬ 
ing ohmic losses. We have shown in a previous work [8j 
that this contribution is not taken into account when 
using the so-called plasma prescription in the Lifshitz- 
Matsubara sum formula which has then to be corrected 
accordingly. However, contrary to non-magnetic materi¬ 
als, the inclusion of a frequency-dependent permeability 
p,(w) can lead, in principle, to either a repulsive or an 
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attractive contribution from those Foucault modes. 

In summary, we have extended the careful analysis 
started in [Sj to magnetic materials. We have been able 
to prove that there is no discontinuity between the Drude 
model, which describes the low frequency behavior of 
metals satisfactorily, and the plasma model, which does 
not, provided the latter is understood as the limit of 


the former where the dissipation is taken to zero. Un¬ 
fortunately, most experiments seem to favor the non- 
dissipative prescription in which the dissipation is ex¬ 
cluded altogether. This problem, called the “Casimir 
puzzle” in the introduction is seen in experiments involv¬ 
ing magnetic materials as well as non magnetic ones and 
it still remains to be solved. 
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